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1. Introduction 

The accurate measurements [jl], |2| of the total cross section near threshold of the 
process 

p+p^p + p + ir° (1) 

have stimulated many theoretical investigations g, f§ §, §, [7|, |, |, [U], 0, 0], but a the- 
oretical explanation of the data is not yet completed. The difficulties in describing the 
process from the low-energy effective field theories such as the heavy baryon chiral pertur- 
bation theory (HBChPT) [|13] are twofold: Firstly, the leading tree order contributions of 
both the impulse (IA) and the meson-exchange (MEC) diagrams are suppressed and the 
contributions from sub-leading IA and MEC diagrams are almost canceled off. Secondly, 
the momentum of the process at threshold is of the scale of ~ ^m n mN where and mjy 
are the pion and proton mass, respectively, which is considerably bigger than the usual 
characteristic scale, m n . Thus contributions from higher order operators - which in general 
are controlled only poorly - can be of non-negligible importance calling for the calculation 
of higher order contributions in the expansion. These are the two principal reasons for 
questioning the predictive power of the HBChPT approach (in its standard form of the 
Weinberg's counting ruleQ]) to the reaction in question. 

According to the chiral filter mechanism [15, 16 1, the processes such as isovector Ml 
transitions and axial-charge weak transitions, which are dominated by one-soft-pion ex- 
change terms (i.e., current algebra), are chiral-filter-protected so that corrections to the 
leading order terms are suppressed and can be systematically controlled by chiral pertur- 
bation theory. When a process is not dominated by soft-pions for reasons of symmetry 
and/or kinematics, then it is unprotected by the chiral filter and consequently higher or- 
der terms (involving short-range ones) become non- negligible, making a systematic chiral 
expansion difficult, if not impossible. The suppression of the leading order contribution 
and the substantially cancellation between the sub-leading contributions of the process @) 
make this process quite similar to other chiral-filter unprotected cases such as the isoscalar 
Ml and E2 matrix elements in the polarized np-capture process, n + p — > d + 7, discussed 



in |lq| and the solar "hep" process studied in [17|. However unlike these processes which 
can be calculated with some confidence because of the small momentum involved, the 7r° 
production process which involves a relatively large momentum has the second (kinematic) 
condition which makes the calculation even more difficult. We hope in this paper to shed 
some lights on these issues which have not been fully explored up to date. 

HBChPT is a consistent and systematic low-energy effective field theory, whose ex- 
pansion parameter is Q/A x , where Q is the typical momentum scale involved and/or pion 
mass, while A x ~ ~ 4irf n is the chiral scale, f T ~ 93 MeV the pion decay con- 
stant. Due to the large momentum scale of the process, Q ~ ^/m^rriN, the convergence 
of the chiral expansion |, I and non-relativistic treatment flTT| , [l2| of the process have 
been questioned. However, there have been many attempts to apply the theory because 
the scale Q is still smaller than A x , (e.g., Q/A x ~ ^Jm^/mN — 0.4). Some tree-order 
calculations g, || [5], |8| and partial one-loop calculations || have been reported. Quite 
recently, Dmitrasinovic et al. |l^] analyzed all relevant transition operators of one-loop 
order and found that important contributions may rise from the one-loop diagrams. In 
this paper, we report the result of our calculation of the reaction cross section by using 
"hybrid approach" |l8| with those transition operators including the initial and final state 
interactions (ISI and FSI) , and short-ranged contributions whose coefficients are deter- 
mined through the resonance-saturation assumption [pTsfj . Indeed, we confirm here that the 
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one-loop contribution is actually quite important, which is not surprising for those pro- 
cesses that are chiral-filter-unprotected |[(], 17 1. Whether or not a reliable prediction can 
be made in the EFT framework, even though the loop-order corrections are not negligible - 
as in the case of the np-capture and the "hep" process - will be the subject of our discussion. 

2. Hybrid approach for the pion production 
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Figure 1: Generic diagram of the pp — > ppir . The empty blobs represent initial- and final- 
state interactions (ISI and FSI), and the shaded blob the two-nucleon irreducible transition 
operators. The momenta in the asymptotic region have bars attached to them. 



The generic diagram for the process is drawn in Fig. |]. In calculating the process, 
we adopt the so-called "hybrid approach", where the two-nucleon irreducible transition 
operators (shaded blob in the figure) are obtained within HBChPT, and all the reducible 
parts are embedded into the phenomenological wavefunctions. The justification for this 



procedure as a viable EFT was given flq , [17 , |20f| and will not be expounded here again. We 
note that it is consistent with the spirit of Weinberg's original proposal of applying EFT to 
nuclear physics |14j . It is well known that, near threshold, the reaction is dominated by the 
transition between | 3 -Po) (initial) and I 1 So) (final) states. With other partial contributions 
neglected, the shape of the total cross section is well reproduced by taking into account 
only the phase factor and the final state interaction O. So we calculate the transition 
amplitude at the threshold kinematics, = (m^, 0) for the four momentum of the emitted 
7T°. According to Sato et al. 0], this approximation reduces the cross sections by up to 10 
% in the range < \q\ < 0.4 m n . 

For the discussion on the kinematics, we denote the incoming (outgoing) momentum 
of the j-th proton by p^ (p'j)- in terms of the relative momentum p and p' ', they can be 

written as p\ = p, P2 = ~P, p'\ = p' an d p'2 = —p' m the CM frame. For convenience, we 
also define the momenta transferred, kj = {pj — p'j)* 1 . The momentum conservation then 
reads 

K + = = (m*, 0). (2) 

Hereafter we will put bars on the momenta defined in the asymptotic region where particles 
are on-shell^; \kj\ = \p\ = ^Jm^m^ + m£/4, \p'\ = and kj = for j = (1, 2). 

The transition operator depends on the energy-transfer, kj = p'j—p'j , which brings up 
an "off-shell" ambiguity [||, 0|. A common method (used in this work as well) is to assume 
the "fixed kinematics approximation" (FKA), which replaces kj by the on-shell energy 

#1 The on-shell condition in HBChPT is given by p 2 j jTe \ — m 2 N — p| + 2mjvu • pj =0, where Pj rcl is the 
usual (relativistic) four- momentum: Pj Icl = mjvt; M + Pj where u M is a constant vector with v 2 = 1. 
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transfer kj = On the other hand, it has been pointed out that different treatments of 
the off-shell behavior can produce significant differences. Sato et al. M suggested to use, 
instead of FKA, the equation of motion approximation (EMA), replacing p° by p 1 /2m at and 

similarly for p'°; thus k® EMA = k® EMA = p 2 ~^ . With this assumption, they carried out a 
momentum-space calculation and obtained quite different results from those obtained with 
FKA. While the EMA can be regarded as an approach containing more dynamics than the 
FKA, neither is free from drawbacks. For instance, it is to be noted that the EMA violates 
the energy conservation, Eq. @. 
3. The total cross section up to one-loop order 

Weinberg's counting rule |L4| dictates that a Feynman diagram with v = 2(1 — C) + 
2L + J2i H is °f order of (Q/A X Y where C (L) is the number of the separate pieces (loops) 
and the index i runs for all vertices. The vi denotes the chiral order of the i-th vertex, 
v = d + n/2 — 2, where n is the number of nucleon lines and d the number of derivatives 
or powers of m n involved in a vertex. As the momentum scale encountered in the present 
reaction is y / m^m/v, a separate counting on (m 7r /A x ) and {y/m^niN j 1 A x ) can be envisaged 
as the order counting for better convergence, as proposed by Cohen et al. |], |2lfl . We 
continue to use Weinberg's counting rule because m^/ A x is still smaller than one. 

Here we wish to remark about the nucleon propagator. Since the nucleon kinetic term, 
p 2 /(2mjv), is of order of m„-, one may think that it should be included in the leading-order 
nucleon propagator, l/(p° + ie). However, this is not the case for irreducible loops^j. The 
energy component of the loop momentum p° can pick up a pole of the pion propagator, which 
has the large momentum scale, and leads to p° ~ ^/m^m^. As a result, p° » p 2 /(2m^), 
and therefore the kinetic energy term can be treated perturbatively. 
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Figure 2: Diagrams for the pion production up to one-loop order. The cross diagrams of 
(b5), (b6), (dl) and (d2) are omitted. 

Feynman graphs relevant to the process are drawn in Fig. where vertices of v = 1 

* For reducible diagrams, we havep ~f/{2m N ). This means that we should include the kinetic term in 
treating the reducible loops, which we always do in solving Schrodinger equations. In other words, the fact 
that p° ~ p°/(2mjv) in reducible diagrams is a general property of low-energy two- and many-body systems 
(one can recall the virial theorem here) , not a particular consequence of the large momentum scale. 
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are marked by "X", and those of u = 2 by a filled circle. Due to the fact that, since 
the transition operator in between the initial state interaction (ISI) and the final state 
interaction (FSI) is of off-shell, a loop integral can have an off-shell singularity that breaks 
the symmetry p2| . We therefore employ the background field method (BFM) [23] which 



preserves symmetries even for off-shell quantities such as the Green's function and effective 
potential. Among the 14 diagrams, eight graphs (a, bl, b2, b3, b6, c, dl, d2) contribute 
here: Graphs (b4), (d3), (d4), (d5) vanish at threshold, (b5) is purely isovector- vector in the 
sea-gull vertex, and the contribution from (d6) is identically zero due to isospin symmetry. 
Note that (b3), (b4), (b5), (d3), (d6) depends on the representation of pion field, but their 
sum does not [pif| . We have not considered the recoil diagrams (Fig. 4 in Ref.Q), which 
are reducible in our scheme and to be absorbed in the wavefunctions. It has been reported 
that the recoil contribution is small Hl-F^ 



The HBChPT Lagrangian is expanded as 

C = Y J ^u = C + C l +C 2 + ---. (3) 
The Lagrangian relevant to our calculation reads 

C = N[iv-D + 2ig A S-A}N + f 2 Tr(iAHA^ + ^\, (4) 



£i = N 



-D 2 + g^^ ■ A, S ■ D} + Cl Tr( X+ ) + ( - 4c 2 J (v • A) 2 - 4c 3 A • A 



2m N m N L 1 ' 1 ™ \2m N 



N. 

(5) 



C 2 = N— — \-^-v AA ■ D - Ac 2 Tr(v ■ AA^)D. ) N + h.c. 
2m. jy \2mN J 

+ Ta \d { ?N[v • AS • D - S- D v ■ A]NNN 

+d { pNv ■ AS^NN[D^- D^N + df N[v • AD^- v • A]NNS^N 

+d ( i ) Nv ■ ANN[S -D-S-D]N + d^ie^v^Nv ■ AS^D p NNS v N 

+ND p v AS^NNSyN - Nv ■ AS^NN S y DpN - Nv ■ AS^NN Dp S„N] } , (6) 

where, in the absence of external fields, = + T^, = ^[£t,c^£], A M = ,9^}, 

X+ = + £x^£> X = m \ an d £ = exp (i^jf-j, V* = (1,0) is the four- velocity vector, 

= (0, ^) is the spin operator and g A is the axial-vector coupling constant, c, and d\ 

are low energy constants which cannot be fixed by the symmetry. The values of a- are 
fixed by using resonance saturation below. For the values of the c's we use those obtained 



#3 In time-ordered perturbation theory (TOPT) the recoil diagrams appear as irreducible. The leading 
order recoil contribution should be, however, removed provided that we are using the usual static one-pion- 
exchange (OPE) potential instead of the non-local one of the TOPT, as Weinberg observed in connection 
to the three-nucleon potential JliJ]. The recoil contribution given in Ref. || is one higher order than the 
mentioned leading order one, and can be obtained by the reducible OPE diagrams (Fig. 2(a,b) of Ref. ||) 
with Feynman pion propagators, subtracted by the same diagrams with static pion propagators. Imposing 
the asymptotic condition for the external lines of the diagram, however, one can show that the recoil 
diagrams give us vanishing contribution at this order. It then implies that the recoil contribution given in 
U is reducible in TOPT as well as in our covariant scheme. Thus including the contribution causes a double 
counting in both schemes, provided that all the reducible diagrams are embodied in the wavefunctions. 
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by Bernard et al. 



ci = -0.93 ±0.10, c 2 = 3.34 ± 0.20, c 3 = -5.29 ± 0.25 [GeV^ 1 ]. 



(7) 



Note that all the degrees of freedom other than nucleons and pions have been integrated out 
from the Lagrangian, their roles encoded in the coefficients of the Lagrangian. In fact, there 
has been a claim that the direct treatment of A(1232) is important Q, and the two-body 
A-contribution substantially cancels the leading-order one-body contribution. The same 
cancellation has been observed in []|], where A(1232) is accounted for only in the constants 
C2 and C3 [26 1, so that we may consider that the role of the A(1232) is reasonably reproduced 
at least at tree order. However, this may not be the case in the one-loop order, since the 
energy scale of loops p° ~ yjirn^m^ is about the same size of the delta-nucleon mass gap. 
Therefore, to be more satisfactory, the A(1232) should be included explicitly; we leave this 
extension for future work. 

The most general form of the transition operator effective at threshold reads 



O = Crf + t1) - B 2 ) ■ PA {1) + kA {2) P s =i 



(8) 



where P = p + p' and k = p — p', and r,- (Bj) is the isospin (spin) operator of the j-th 
nucleon, and Ps=i = \{B\ • B 2 + 3) the spin-1 projection operator. In deriving the above 
equation, we have used B\ x B 2 = i{B\ — B 2 )P° , where P a = ^(1 + B\ ■ B 2 ) is the exchange 
operator in spin space. 

The An\ receives contributions solely from the impulse diagram (a), 



A (i) 



gAm n 



0) 



8m N U 

Note that all the higher order corrections at threshold are already included in the above 
equation when the physical values of gA and other parameters are used. On the other hand, 
the A( 2 ) receives contributions only from two-body graphs. The contributions from the 
one-pion-exchange (OPE) diagrams including the vertex loop corrections can be written as 

1 



A {2) 



9A 

4/tt k 2 



mi 



-T nN (k,P), 



(10) 



where k^ = k± = (k°, k). Here and hereafter, we use the FKA; k° = k® = k 2 = m 7r /2. Up 
to one-loop accuracy, the T n N is given as 




1 / 

2ci + - I c 2 + c 3 



32mjv 



-7T+ C2 



9\ 



16mjv 



k ■ P 



m" 



32x/3tt/4 



327T 2 /^ ^0 



x(l — x)k 2 + 3D(x, 0, q — k) n(x,0,q — k), 



(11) 



with 



n(x, uj, k) 
D(x, uj, k) 



x 



k° 



y/D(x,LJ, k) 



7T 

arctan- , 

K 2 y/D(x,w,k)_ 





(12) 

(13) 

#4 The values of Cj's have been updated in |^]. We have checked that our results do not depend on these 
parameter set much, less than ~ 10 % in total cross section. 



x(l — x)k — (xk — uj) 
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The four terms in Eq. (11) correspond to (bl), (b2), (b6), (b3) diagrams in Fig. [2|, 
respectively. Among the two-pion-exchange (TPE) graphs drawn in Fig. 2{dl — d6), only 
(dl) and (d2) contribute to Am 



(27T 



9Am n 



® 128^/5 Jo 



1 — X X n 

dx i — n_ + ~n+ + 3^ 



1.2 

(9- 10x)ln[l - x(l - ar) — =■] 

mz 



^^Qlj ( 1 3^ J /c r / \ i 

+4(5x — 6)ra_ + 12no — 4(5x — 3)n + — [(1 — x)n_ — no + xn + \ 



mi 



,(14) 



where n± = n(x, ±m n /2, —k), uq = n(x, 0, —k), and similarly for the functions with a bar, 
and 



n(x, u, k) 



D(x, u, k) 



mt 



n(x, uj, k). 



(15) 



Finally we have the counter term contributions, drawn in Fig. 2(c), 



with 



A CT _ 9Am n ( 2 ) 
« " 128vr2/r 1 



4 2) +4 2) 



,(2)x ,CT _ J3Am^,(2)R 



(16) 



7 (2) 



1 /^tt^ 
7 + In 



(17) 



(2) F? - 

where dg is the renormalized constant in (4 — 2e) dimension. 7 = 0.5772 • • • and p is the 
renormalization scale of the dimensional regularization. In coordinate space, the counter 
term contributions are accompanied by the delta function. While the initial pp state is in 
P-wave, these zero-ranged operators are effective due to the derivative operators residing 
in Eq. (||). In principle, the d^ 2 ^s should be determined from experiments, a task which 
is however not feasible due to the lack of available data. We instead fix them from the 
resonance-saturation assumption, taking into consideration the omega and sigma meson 
exchanges || 



A 



CT 
(1) 



9Am n / g% 



+ 



.91 



T m 



N 



mf- 



ACT 
A {2) 



9Am n g£ 



4/ 7r m : 



N 



11 r 



(18) 



where g a (g w ) and m a (m u ) are the scalar (vector) coupling and the mass of the sigma 
(omega) meson, respectively. Using the values |J7| g a = 10.5, m a = 508 MeV, g u = 10.1, 



and m w = 783 MeV, we have df ] - d% ] + df ] - df ] = -7.70 and df )R 
we have imposed the zero- momentum subtraction scheme to the A^y 
The total cross section is given by 



1 2?F / A IT|2 

a = / dp\T\ 

4wiab 



with dp 



(2tt)< 



-rn 



>\o\ 2 Pfd\q\, 



-4.78. Note that 



(19) 



where i>i a b is the incident velocity in the lab frame, T the transition amplitude, and dp the 
phase factor; m p is the proton mass, pf the magnitude of the relative three-momentum of 

* 5 One can easily distinguish the contributions from (dl) and (d2) diagrams by noting that the contributions 
of (dl) and (d2) are proportional to qa and g A , respectively. 



7 




Figure 3: Total cross sections up to the one-loop order in HBChPT. The wavefunctions 
are from RSC (solid curve) and HJ (dashed curve) potentials. Experimental data are also 
plotted. 



the final pp state, and q the momentum carried by the outgoing pion in the CM frame. The 
upper limit of \q\, \q\ ma , x , is set by the initial total energy. 

In calculating the amplitude T, we use the wave functions obtained by solving 
Schrodinger equation with phenomenological potentials. In this work we use the Reid soft 
core (RSC) and Hamada-Johnston (HJ) potentials [28, 29|.^ We write the initial and final 
state wave functions as 



*<(r) = iV2V^^e l5 i\ 3 P ), *,(r) = v^^e i5( f 



So), 



(20) 



with the normalization condition as iix(r) r z!^ ) i sin(pr — ttL/2 + 6l). The amplitude T is 
given as 



16vr 



'm n Jo 



dr 



A(iQ {uqu^ - uiu' H — j - A(2)-^ (uom) 



(21) 



where ul = Ul(t), u' l = ^rUi(r), and Aq) = Ay)(r) are the transition operators in 
coordinate space, 



d 3 k 

(>)3 



Ar-k 



A {j) (k). 



(22) 



4. Numerical result and Discussion 

In Fig. 3 we plot the total cross section vs. n = \q\ ma , x /m n . Our theoretical predic- 
tions with loop corrections come close to the experimental data. This clearly points out the 



* 6 In H it has been reported a large sensitivity of the cross section to the short-range part of the potential. 
In our work we will study this feature by introducing a cutoff. 
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Figure 4: Total cross section of our result for various cutoff r c with the wave-function from 
RSC potential. The contributions from the zero-ranged operators are also included with 
replacing 5{r) by S(r — r c ). 



importance of the one-loop contribution. This aspect has been noted before while mention- 
ing chiral-filter-unprotected cases. Note that since the HJ potential has a hard-core whose 
radius is about 0.5 fm, the zero-ranged counter-term contributions are identically zero. 
This corresponds to a "hard-core regularization" of the short-distance physics encoded in 
the counter terms. This was referred to in |l6| as "HCCS". For the RSC potential, the 
counter-term contributions are not zero. To properly calculate the counter-term matrix 
elements in a consistent regularization scheme, one has to resort to what was referred to in 
[16] as "MHCCS" which amounts to replacing the delta function and the two-body part by 

6{r) ^ § (r - r c ), ljj 27r (r) -> 6{r - r c )A^(r), (23) 

The corresponding results with the RSC wavefunctions are given in Fig. 4. As can be 
seen in Fig. 4, the results are mildly sensitive to r c . This mild sensitivity - which may be 
taken as a justification of the short-distance regularization procedure used here - was also 
observed in [16, 17] in different contexts. 

A serious question here is the convergence of the perturbation series. It is difficult to 
see whether or not the series actually converges without computing higher order terms. But 
going to higher order is not feasible at the moment. To gain a rough idea from the present 
calculation, we may consider the break-down of the transition amplitude into individual 
contributions. For convenience of comparison, we measure correction terms relative to 
the IA amplitude calculated with the RSC wavefunctions at |(f|max = 0.1 m n . Then the 
leading MEC (coming from OPE) is —0.8 in units of the IA term, leading to a substantial 
cancellation. The one-loop order MEC contribution is (—1.0) + (—1.2) = (—2.2), where 
the first (second) term represents the individual contribution from OPE (TPE) diagrams. 
Summing up the whole (including IA) contribution, we get —2.0, that is, twice the IA 
contribution with the opposite sign. This shows that there is no convergence to the order 
considered. This feature is shared by other chiral-filter-unprotected cases; in both the 
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polarized np-capture and the "hep" process, a similar pattern is seen. There is however 
one important difference. In the latter cases, the coefficients of the counter terms can be 
reduced to one effective constant and, when this constant is fixed by experiments, can be 
more or less accounted for higher order terms that are not computed (such as the magnetic 
moment of the deuteron in the case of polarized np-capture and the triton beta decay in 
the case of "hep"). This amounts to terminating the series by fiat at the order considered 
by the counter terms. One can think of this as choosing - as one does in gauge theory 
calculations - an optimal regularization scheme so as to account for higher-order terms that 
are not computed explicitly. This freedom of choosing the optimal regularization scheme is 
not available in the tt° production case since the counter term is entirely fixed by means of 
the resonance saturation and hence cannot "mock up" higher-order terms which may not 
be negligible. 

Next, we discuss the effect of ISI and FSI in studying the convergence of the 1/mjv 



correction of the C2 term. In [12] the 1/mjv contribution was reported to become roughly 
twice as large as the leading C2 term; this result was obtained by employing a simple 
"factorization ansatz" which dictates that the nucleon legs of the production operator are 



imposed to be on mass-shell. This is confirmed in [ 10 1 : The l/mjv correction, (b2) diagram 
in Fig. ^, is about 260 % of the leading MEC contribution (bl) when the ISI and FSI are 
ignored. We would obtain the same result if we exclude the ISI and FSI, that is, imposing 
the on-shell condition k ■ P = m^nriN + in Eq. ([TT]) . However, after taking account of 

the ISI and FSI we find the matrix element of (b2) is only about 20 % and 40 % of the (bl) 
evaluated with RSC and HJ potential, respectively. Consequently, the 1/mjy correction 
does not show any anomalous behavior in our calculation. The similar pattern of softening 
the l/m jv correction of the Weinberg- Tomozawa term due to the ISI and FSI was observed 



in [30 1 as well. 

We want to remark on Sato et aVs observation [|3~l| [: In a momentum space calcu- 
lation with the half-off-shell wavefunctions, they found significant contributions from the 
momentum region even larger than the chiral scale A x ~ 1 GeV for some Feynman di- 
agrams. Regarding this apparent disagreement between their cutoff-dependence and our 
r c -independence, we should note that the observed r c -independence in our theory is in the 
total amplitude, not in individual terms. In fact if we evaluate the contribution for each 
diagram, a considerable r c -dependence (though not as radical as theirs) is observed in our 
calculation also. 

In summary, we confirm quantitatively that the two pion exchange one-loop contri- 
bution is quite important. We also find that there is no strong dependence on the short 
range interaction. This result suggests that massive degrees of freedom (not included ex- 
plicitly in our study) may have been properly incorporated into the diagrams and counter 
terms investigated in our study. The apparently poor HBChPT convergence however may 
indicate that significant contributions from higher-order operators may be left unaccounted 
for in our scheme. Furthermore, the validity of the FKA should be examined. Clearly more 
work is required including comparison with the other counting rules ||, studies of next 
order terms, and the verification of the approximations involved. 
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